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Abstract. We discuss how rotation and binary interactions may be related to the diversity of 
type Ibc supernovae and long gamma-ray bursts. After presenting recent evolutionary models 
of massive single and binary stars including rotation, the Tayler-Spruit dynamo and binary 
interactions, we argue that the nature of SNe Ibc progenitors from binary systems may not 
significantly differ from that of single star progenitors in terms of rotation, and that most long 
GRB progenitors may be produced via the quasi-chemically homogeneous evolution at sub-solar 
metallicity. We also briefly discuss the possible role of magnetic fields generated in the convective 
core of a massive star for the transport of angular momentum, which is potentially important 
for future stellar evolution models of supernova and GRB progenitors. 



1. Introduction 

Rotation influ ences not only the evolut i on of massive stars, but also the ir supernova 
(SN) explosions (jMaeder fc Mevnet llioool : iHeger. Langer fc Wooslev llioooh . In particu- 
lar, recently many asymmetric supernovae with unusually large energy (broad-lined SNe 
or hypernovae) have been discovered (e.g. iMazzali et al. 1120071 ). which shows evidence 
for rapidly rotating progenitors. The most spectacular example may be long gamma-ray 
bursts (GRBs), which are generally believed to be produced by deaths of some massive 
stars retaining extremely large angular mom enta in their cores (7 1 



16 cm2 s 



Wooslev 1ll993t MacFadven fc Wooslev 1999t see, however. iDessart et al. I [20081 ). Inter- 
estingly, such energetic core-collapse events seem to only occur in Wolf-Rayet (WR) stars: 
all of the broad-lined SNe/Hypernovae and the supe rnovae associated with lo ng GRBs 
have been observationally identified as Type Ic (see Wooslev fc Bloonil |2006| for a re- 
view). This raises the question which WR stars can produce broad- lined SNe Ic or long 
GRBs while most WR stars die as normal SNe Ibc. Here we present recent evolutionary 
models of massive stars that include binary interactions and the transport of chemical 
species and angular momentum via rotationally induced hydrodynamic instabilities and 
magnetic torques, and discuss how rotation and binary interactions may be related to 
the diversity of SNe Ibc and long GRBs. 



2. Single star models 

Redistribution of angular momentum and chemical species in a rotating st ar occurs by 
rotationally induced hydrodynamic instabilities and magnetic torques (See iTalon 2007 
for a review). Eddington-Sweet circulations, shear instability and Goldreich-Schubert- 
Fricke instability among others have been considered in previous non-magnetic models 
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Figure 1. Mean specific angular momentum of the star and tlie innermost 1.4 M© and 3.0 M© 
as a function of tlie evolutionary time for 18 Mq (left panel) and 60 Mq (right panel) models. 
The time span for the helium core contraction is marked by the color shade as indicated by the 
label. 



(iMaeder fc Mevnet |[2000l : iHeger. Langer fc Wooslevl '2000': 'Hirsc hi. Mevnet fc Maeder 



20041 ) ■ and the so-called Tayl er-Spruit dynamo dSpruit 2002.) has been implemented in 
recent magnetic mode ls (e.g. iHeger. Wooslev fc Spruit II2005I : Maeder fc Mevnet 2005 



Yoon fc Langer 1 1200,^ . In non-magnetic models, it is shown that the buoyancy due to 



the chemical gradient at the interface between the core and the envelope largely prohibits 
the considered rotationally induced hydrodynamic instabilities from transporting angular 
momentum. The amount of angular momentum retained in t he core at the pre-supernova 



stage is thus close to its initial value eve n at solar metallicity (|Heger. Langer fc Woosley 



20001 : iHirschi. Mevnet fc Maeder Il2004( l. Most massive stars are predicted to die with an 
enough amount of angular momentum in the cores to produce long GRBs via formation 
of millisecond magnetars or coUapsar (i.e., j >~ 10^^ cm^ s~^), given that a large frac- 
tion of young massive stars in our Galaxy and Small/Large Magellanic Clouds are rapid 
rotators (e.g. lMaeder fc Meynet1l2000HMokiem et aLlliooalHunter et al. Ibooj ). On the 
other hand, in magnetic models adopting the Tayler-Spruit dynamo the core is effectively 
spun down by magnetic torques, and the predicted spin rates of white dwarfs and young 
neutron stars are smaller by two orders of magnit ude than those from non-magnetic mod 



els, w hich can better explain observations (Heger. Wooslev fc Spruit 1 120051 : ISuiis et al. 



20081) . 



Fig. [T] shows the evolution of the core angular momentum in the magnetic model 
sequences with Minit = 18 M0 & Urot,init = 144 km s~^, and 60 M© & Wiot,init = 
186 km s^^. In the sequence with Minit = 18 M©, the core loses a significant amount of an- 
gular momentum during the helium core contraction phase where a strong degree of differ- 
ential rotation between the core and the envelope appears. A similar effect is also observed 
during the CO core contraction phase. In the sequence with Minit = 60 M©, spinning- 
down of the core during He-core contraction becomes less significant as the star loses the 
hydrogen envelope, which leads to a smaller moment of inertia of the envelope. However, 
the core is further spun down by loss of mass due to LBV/WR winds during core He burn- 
ing. At the neon burning stage, both stars retain a similar amo unt of angular momentum 
in the ir cores as shown in Fig. [1] In fact, the calculations by lYoon. Langer fc Norman 



dlooi) show that magnetic models give < ji.4 >~ 2. ..3 x 10^"* cm^ s ^ for different initial 



metallicities, masses and rotational velocities, implying that most massive stars including 
Type Ibc progenitors should die with a similar amount of core angular momentum. This 
conclusion remains the same even for binary stars, as discussed below. 

An exception is the case for the so-called chemically homogeneous evolution. If the ini- 
tial rotational velocity is exceptionally high and if metallicity is sufficiently low, chemical 
mixing by Eddington-Sweet circulations may occur on a time scale even smaller than 
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Figure 2. Final fate of our rotating massive star models at four different metallicities {Z — 
0.004, 0.002, 0.001, & 0.00001), in the plane of initial mass and initial fraction of the Keplerian 
value of the equatorial rotational velocity. The solid line divides the plane into two parts, where 
stars evolve quasi-chemically homogeneous above the line, while they evolve into the classical 
core-envelope structure below the line. The dotted-dashed lines bracket the region of quasi-ho- 
mogeneous evolution where the core mass, core spin and stellar radius are compatible with the 
coUapsar model for GRB production (absent at Z=0.004). To both sides of the GRB production 
region for Z — 0.002 and 0.001, black holes are expected to form inside WR stars, but the core 
spin is insufhcient to allow GRB production. For Z = 0.00001, the pair-instability might occur 
to the right side of the GRB production region, although the rapid rotation may shift the pair 
instability region to larger masses. The dashed line in the region of non-homogeneous evolution 
separates Type II supernovae (SN II; left) and black hole (BH; right) formation, where the min- 
imum mass for BH formation is simply assumed to be 30 M©. From lYoon. Langer fc Norman I 
(1200^) . 



the nuclear time scale. Quasi-homogeneity of the chemical composition of a star is thus 
ensured on the main sequence and the star is gradually transformed into a massive WR 
star, avoiding the giant phase that would result in strong braking down of the core. All 
of the necessary conditions for producing long GRBs - massive core to make a black 
hole, removal of hydrogen envelope and retention of a large amo unt of angular momen- 
t um in the core - thus c an be fulfilled by this type of evolution ( Yoon fc Langer I [2005 : 



Wooslev fc Heger II20Q6I ). This chemically homogeneous evolution scenario (CHES) fa- 



vors low metallicity environment for producing long GRBs (Fig. [2]) and predicts a higher 
ratio of GRB to SN rate at highe r re dshift, which should be tested by future observations 
(|Yoon. Langer fc NormaTlbOOd cf. lKistler et al. Il2007f ). 



3. Binary star models 

A significant fraction of SNc Ibc may be produced in close binary systems (e.g. 
Podsiadlowski. Joss fc Hsu 1992). An example is given in Fig.[3]that shows the evolution 
of of the primary star in a close binary system with Pjnit — 4 days, Mpi.imary,init — 18 Mq, 
and Msccondary,init = 17 Mq . Both stars are tidally synchronized early on the main se- 
quence. Once the primary star fills the Roche-lobe radius, it loses about 7 Mq during 
the Case A mass transfer phase, and additional 3.5 Mq later during the Case AB mass 
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Figure 3. Left panel Evolution of the internal structure of the primary star in a binary system of 
Mprimary.init = 18 Mq , Mgecondary ,init = 17 Mq and Pinit = 4 days, from zero age main sequence 
to neon burning. Right panel Mean specific angular momentum of the innermost 1.4 Mq and 
3.0 M0 of the primary star considered in the left panel as a function of time. The time span 
for Case A or Case AB mass transfer phase is marked by the color shades as indicated by the 
labels. 



transfer phase, becoming a 4 Mq WR star. The core loses angular momentum mostly 
during these Case A and AB mass transfer phases as shown in Fig. [31 The amount of 
angular momentum in the core at the neon burning stage turns out to be very similar 
to those in single star models (see Fig. [Ij. In fact, we find that mean specific angu- 
lar momentum of the innermost 1.4 at the final evolutionary stage of a primary 
star in a binary system does not change much according to different initial parameters 
(primary mass, mass ratio and orbital separation): it remains within a narrow range of 
2. ..3.5 X 10"'^'' cm^ s~^ regardless of the detailed history of binary interactions, as long 
as the tidal synchronization is not unusually strong (Yoon, Woolsey & Langer 2008, in 
prep.). This implies that the nature of most binary star progenitors of SNe Ibc may not 
much differ from that of single star progenitors, in terms of rotation. 

The evolution of secondary stars has not yet been well understood. In particular, 
it sensitively depends on the uncertain e fficiency of semi-convec tion whether the mass 
accreting star may be rejuvenated or not ( Braun fc Langer"lll995l ). If a rather large semi- 
convection parameter is adopted, rejuvenation can significantly weaken the chemical 
gradient between the hydrogen burning core and the envelope, in favor of rotationally 
induced chemical mixing. As the secondary is spun up to the critical rotation by mass 
accretion, even the chemically homogeneous evolution can be occasionally induced if 
metallicity is sufficiently low and if the secondary is not strongly spun d own by the tidal 



synchronization after the mass accretion phase ( Cantiello et al. 2007t ). The secondary 



will eventually die as a GRB after traveling from a few to several hundreds PCs away, 
if the binary system is unbound due to the supernova kick as a result of the explosion 
of the primary. This scenario may ex plain the recent obser vational evidence that some 
GRBs are produced in runaway stars (jHammer et al. hood ). 



Other types of binary interactions may also lead to formation of rapidly rotating WR 

stars to produce long GRBs. Tidal spinning-u p of a WR star in a compact binary system 

with a neutron star or a black hole companion (iBrown et. al. l2000l: l llzzard. Ramirez-Ruiz fc Tout 



2004l:lvan Putten. M.H.P.M. Il2004l : lvan den H euvel fc Yoon I I2007I ) and merger of two he 



Hum cores in a common evenlope (jFrver fc Hc gcr 2005.) have been recently suggested 

among others. It remains uncertain, however, that such binary sys tems could explain the 

observed GRB rate. For example, recent stellar evolution models bv lDetmers. Langer fc Podsiadlowski 



(2008 ) show that a merger of the WR star with the compact object, which is not sup- 
posed to produce a classical long GRB, is the most likely outcome in the former case. 
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Figure 4. Mean radial fields Bi (r, 9) on the meridional plane in a 12 Mp, rotati ng star on the 
main sequence in a MHD simulation with a 3-D anelastic code (jGlatzmaier II 19841 ) . The adopted 



angular velocity is 10 Rad s 
units) is the convective core. 



The inner region of r ^ 7 x 10 cm (r 0.5 in the code 



Further detailed evolutionary model s are certainly needed for observationally testing 
different evolutionary scenarios (e.g. Ivan Marie et al. On the other hand, it is 

puzzling why no GRB-associated SNe lb h ave been observed yet w hile most GRB pro- 
genitor scenarios predict their existence (e.g. lYoon. Langer fc Norm an 2006). Within the 
C HES, this puzzle might be solved if one considered anisotropic mass loss, as discussed 
in lMevnet fc Maeder I (|2007[ ). 



4. Concluding remarks 

Our stellar evolution models including the Tayler-Spruit dynamo indicate that most 
SNe Ibc progenitors should explode with a similar amount of angular momentum in their 
cores, regardless of their single or binary star origin. This is due to the self-regulationary 
nature of the Tayler-Spruit dynamo. Loss of the hydrogen envelope due to stellar winds or 
mass transfer results in both removal of angular momentum from the core and weakening 
of the core braking by the extended hydrogen enveloped due to magnetic torques, and 
vice versa. Therefore, most different pre-supernova evolutionary paths may not contribute 
much to the diversity of Type Ibc supernovae in terms of core angular momentum, 
although diff erent iron core masses, and thu s different spin rates of young neutron stars 
may result ((Heger. Woosley fc Spruit GRB progenitors, which require unusually 

large angular momenta, may undergo the chemically homogeneous evolution, which may 
not be unusual at low metallicity. On the other hand, recent obs ervations indica te that 
not all broad- lined SNe Ic are associated with long GRBs (e.g. ^Modiaz et al.ll2008 ). which 
needs a theoretical explanation in future work. 

Although the predicted spin rates of the stellar remnants from the magnetic models are 
consistent with observations, the validity of the Tayler-Spruit dynamo has been recently 
questioned by several authors ( Denissenkov fc Pinsonneault 12007 : Zahn. Brun fc Mathis 
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20081) . Furthermore, we might have ignored potentially important physical ingredients in 



simulating the evolution of rotating massive stars. These include gravity waves (Townsend 
in this volume), and magnetic fields generated by the convective core . For instance, our 
recent 3-D simulations with an anelastic magnetohydrodynamics code ( Glatzmaier Ill983) 
show that the strength of poloidal fields generated in the convective core in a young mas- 
sive star may amount to several thousand Gauss on average (Fig. and its influence on 
the transport processes might be comparable to what the Tayler-Spruit dynamo predicts. 
This issue will be addressed in Yoon, Woosely & Glatzmaier (2008, in prep.). 

This work was, in part, supported by the DOE Program for Scientific Discovery through 
Advanced Computing and NASA. 
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